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Angular resolved photoelectron spectra of metal clusters have been experimen- 
tally measured for the first time only recently These measurements have been 
performed systematically for sodium clusters in a broad range of cluster sizes. This 
work attracted a lot of attention and was reported practically at all major interna- 
tional cluster conferences because it revealed a very non-trivial behavior of the angu- 
lar anisotropy parameter with respect to photon energy and provided a method for 
probing the angular momentum character of the valence orbitals of free nanoclusters. 
Initial attempts to explain these observations within single particle approximations 
fail completely In this Letter we present a consistent many-body theory for 
the description of angular resolved photoelectron spectra of metal clusters. Jellium 
model formalism is employed. Our calculations demonstrate the dominant role of 
the many-body effects in the formation of angular distributions of photoelectrons 
emitted from sodium clusters and are in a good agreement with experimental data 
reported in The concrete comparison of theory and experiment has been per- 
formed for the photoionization of Na^ and -/Va^g anions being characterized by the 
entirely closed shells of delocalized electrons. 

PACS numbers: 36.40.-c,32.80.Gc,33.60.+q,36.40.Wa, 



The problem of photoionization of metal clusters has been studied intensively during the 
last two decades, see e.g. [2I and references therein. During this period a number of impor- 
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tant achievements have been made. Thus, the plasmon resonances have been observed in 
small metal clusters of nanometer size for different materials, e.g. Na, K, Mg, in fuUerenes 
and in many other cases of study, for references see j^. The structure of these plasmon res- 
onances have been studied in detail. It was established that patterns of plasmon resonances 
in many cases are determined by the deformation parameter of the cluster. This result is 
very important because it provides a method of monitoring the cluster size and its shape. 
A comprehensive review of the results in this field can be found in [s]. 

The detail geometrical structure of mass selected clusters has been also studied by means 
of photoelectron spectroscopy and the density functional theory (DFT). Thus, for sodium 
cluster anions with up to 57 atoms the geometrical structures have been determined with 
high accuracy from the comparison of DFT calculations results with measured photoelectron 
spectra {4, 5]. 

The measurement of photoabsorption and photoionization of small mass selected clusters 
is not a trivial task. Until recently the measurement was fulfilled only for the total pho- 
toionization or photoabsorption cross sections of various cluster targets. In the recent work 

n 

[11 the first measurement of angular resolved photoelectron spectra of sodium clusters has 
been reported. These experiments have been performed for negatively charged sodium ions 
(anions) in a broad range of cluster sizes 3 < < 147 {l|. These experiments allowed to 
orobe the angular momenta of single electron orbitals in sodium metal clusters. In the work 
11] it was also demonstrated that simple models based on single electron treatment of the 
photoionization process fail to describe the angular anisotropy of photoelectrons emitted in 
the process of photoionization of cluster anions. 

In this Letter we present a consistent many-body theory for the description of angular 
resolved photoelectron spectra of metal clusters. As a case of study sodium cluster anions 
have been chosen. The Hartree-Fock (HF) approximation has been used as a single particle 
theoretical framework. Many-electron correlations have been accounted for within the Ran- 
dom Phase Approximation with Exchange (RPAE). The results of calculations allow one 
to conclude that many-electron correlations play the very essential role in the formation of 
angular distributions of photoelectrons in the process of photoionization of metal clusters. 
This effect is a consequence of the large dynamic polarizability of metal cluster targets, being 
entirely determined by many-electron correlations in the vicinity of the plasmon resonance 
frequencies. 
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In this Letter we consider the most characteristic case of a spherically symmetric cluster 
target and apply the jellium model for the description of the photoionization process. This 
model proved to be well applicable for the description of electronic structure and colli- 
sion processes involving metal clusters and fuUerenes . Cross sections of the collision 

S'ocesses are very sensitive to the correct accounting for many-electron correlations, see 
, and references therein. In this Letter we limit our description by the cluster targets 
possessing the spherical symmetry, which corresponds to the case of the so-called magic 
clusters, i.e. clusters with the entirely closed electronic shells. For the concrete analysis and 
illustration we have chosen the Naj and Na^g magic clusters. 

Let us note here also that the angular anisotropy of photoelectrons emitted in the process 
of photoionization of atomic negative ions has been investigated in sufficient detail both 
theoretically and experimentally (see, for example, review |6|). It was demonstrated that 
the parameter of angular anisotropy /? [3] is very sensitive to accounting for the many- 
electron correlation effects. Thus, it is not a surprise the failure of single particle approaches 
reported in jl| for the case of the photoionization of cluster targets. 

The characteristic features of small metallic clusters, like the shell structure of delocalized 
electrons or plasmon excitations, can be well understood in terms of quantum motion of the 
delocalised valence electrons moving in the field created by themselves and the positively 
charged ionic core Q]. This concept, known as the jellium model for metallic clusters, can 
be also utilized for the description of a metallic cluster anions. 

In this work for the elucidation of the role of many-electron correlations in the formation of 
angular distributions of photoelectrons we have addressed to the most characteristic example. 
Thus, we considered the magic cluster anions possessing the spherical symmetry due to the 
entire closure of their electronic shells. In this case the complicated ionic structure can be 
reduced with the sufficient accuracy to the uniformed spherically symmetric positive charge 
distribution, named below also as ionic core. 

Within the Hartree-Fock (HF) approximation the ground state many-electron wave func- 
tion of the electronic subsystem is defined by the Slater determinant constructed of single- 
electron wave functions fi{x.), being solutions of the system of coupled integral-differential 
equations: 

(T + U^re + Vel)^i{^) = i (l) 

where T is the kinetic energy operator, Ucore is the potential of the ionic core and V^i is 
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the electron-electron interaction including direct and exchange terms. The atomic system 
of units \e\ = rUe = h = 1 is used throughout the paper. 

When calculating the ground state wave functions, the system of equations ([I]) should be 
treated as coupled and to be solved self-consistently. Each equation in ([T]) can be solved by 
an iterative procedure. The eigen-values of these equations are the single-electron energies 
ej. In a spherically symmetric field, the single electron wave functions in ([1]) can be expressed 
as Lpni,ii,mi{^ , where nj,/j,mj are the principal, orbital, magnetic quantum numbers of the 
orbital i. This means that when neglecting the spin-orbital interaction the cluster electronic 
shells with quantum numbers nl consist of 2(2/ -|- 1) electrons. For instance, in the case of 
Ntty and Na^g the electronic configurations of the ground state have the form: Is^ 2p^ and 
Is^ 2p^ 3(i^° 2s^, respectively. 

The wave function of a detached photoelectron can also be obtained from equation 
as a solution with the fixed energy e = Ip + u obeying certain asymptotic behavior 
where Ip is the ionization potential of the initial state and uj is the photon energy. The wave 
functions of excited states have been calculated both in the field of a "frozen" electronic 
cluster core with the created vacancy and with the entirely rearranged residual electronic 
structure of the cluster. The latter case corresponds to the situation of the so-called static 
rearrangement. 

Within the dipole approximation the angular distribution of photoelectrons ionized from 
raZ-state by linearly polarized light can be written in a well-known general form [3]: 

daniiyj) (Tnlioj) 
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[l+(3nl{uj)P2{cOs9)), (2) 



where ct„;(ci;) is the total photoionization cross section from nl state, u is the photon energy, 
P2(cos^) is the Legendre polynomial of the second order, 6 is the electron emission angle 
with respect to the polarization of the incident light. The angular asymmetry parameter /5 
[3] in the single-electron approximation reads as ^ 



6a//(/ + l)di^idi+i cos{Si+i - Si^i)]. (3) 



Here 6i±i is the phase of photoelectron, / is the orbital momentum of the ionized shell 
and di±i is the reduced dipole matrix element, which can be obtained from photoionization 
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amplitude within the dipole approximation {el ± \ \d\nl) by the integration over angles and 
the summation over spin variables. In the length form it reads as 



l> J Pniir)rP^i±i{r)dr, 





(4) 



where Pni{r) and Pei±i{r) are the radial parts of initial and excited states wave functions, 
/> = / + ! for transitions [l^l + 1) and /> = I for the case of {l^l — 1). 

There is a certain threshold value for the angular anisotropy parameter (3 in the cross 
section of photodetachment from negative ions |6| . For / > 1 it is equal to 

/ - 1 



nl 



21 



(5) 



In the case / = the (3 parameter is equal to 2. 



Expression for the angular distribution of photoelectrons within consistent many-body 
theory accounting for many-electron correlations has been derived for the first time in [lO|. 
In this formalism the reduced dipole matrix elements become complex. The expression for 
/3 becomes rather cumbersome and reads as js]: 

1 



/3riz(u;) 



{(^-1)1^, 



+ (/ + 2)|A+i|' + 6V/(/ + l) X 

X [{ReDi^iReDi+i + ImDi^iImDij^i) x 

X cos((5/+i - 5i^i) - 

- {ReDi^iImDi^i - ReDi+iImDi^i) x 

X sin(5i+i - 



/-I 



(6) 



where Di±i are the complex reduced amplitudes obtained from the dipole matrix elements 
{el ± l\D{u)\nl) being solutions of the RPAE equation. In the matrix form the RPAE 
equation reads as j^: 



+ 



{el ± l\b{uj)\nl) = {el ± l\d\nl) + 

^ {t,\D{Lu)\is){t3,el±l\U\t,,nl) 



i4>F 



i3>F 
i4<F 



64 €3 + i5 



(7) 
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where indices and describe the virtual electron-hole excitations, the combined matrix 
element {U) is a sum of the direct and the exchange Coulomb matrix elements (V): 

{is, el ± l|i/|z4, nl) = 
= {i3,el±l\V\H,nl) - {el±l,i3\V\i4,nl) (8) 

In calculations we take into account correlations between all possible transitions nl — >• 
el ± 1. 

Using the formalism presented above we have calculated the angular anisotropy param- 
eter f3 for the photoionization cross section of the magic Na^ and Na^g cluster anions. 
Calculations have been performed in the single particle HF approximation and also with 
accounting for many-electron correlations. 

n 

In Fig. [1] we present comparison of the experimental data from [l[ with the results of 
our calculations of angular anisotropy parameter j3 for the photoionization of 2p shell of 
the Na^ cluster performed in the single particle HF approximation. Different experimental 
dependencies A, B, C in Fig. [1] correspond to the photodetachment of p-electrons from 
different sublevels of the 2p-orbital arising due to its splitting by the crystalline field of 
the cluster. In the paper ll| there was also performed the calculation of the parameter 
/3 in a single-particle approximation with the use of the self-consistent Ekardt potential 
and additional accounting for the image-charge and dipole potentials. Our result here is 
based on Eq. ([3]) with the amplitudes being calculated using the HF wave functions. As it 
follows from Eq. ([5]) the threshold value of f3 for the p-photoelectrons should be equal to 
0. From the results obtained it is clear that within a single-particle formalism it is neither 
impossible to derive the correct threshold value for /5 nor the overall correct behavior of 
(3{uj) corresponding to the experimentally measured dependence. 

In order to elucidate the role of many-electron correlations calculation of the angular 
anisotropy parameter (3 has been performed on the basis of the RPAE equation Eq. ([6]), 
from which the amplitudes Di±i have been derived. In Fig. [2]the results of this calculation of 
(3 for the photoionization of Na^ cluster are presented and compared with the experimental 
results from l]]. Comparison of Figs. [H and [2] demonstrates that accounting for many- 
electron correlations leads to a qualitative change of behavior of the calculated dependence 
and thus the angular distribution of photoelectrons. Fig. [2] demonstrates also that the 
consistent many-body theory leads to the correct threshold value of the parameter (5 and 
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FIG. 1: Angular anisotropy parameter /3 for the photoionization cross section of the Nuj cluster 
versus photon energy; solid line is a result of the HF approach derived in this work; dotted line 

n 

is a result of a single-particle approach used in [l,]; A, B, C are the experimentally measured 
dependencies 
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FIG. 2: Angular anisotropy parameter /3 for the photoionization cross section of the Nuj cluster 
versus photon energy; solid line is a result of the RPAE calculation accounting for many-electron 
correlations; A,B,C are the experimentally measured dependencies from 

improves substantially the agreement of theoretical results with experimental data. 

A similar calculation and analysis of P have been performed for the photodetachment 
from the Na^g cluster anion. In this case we have calculated the wave functions for outgoing 
photoelectrons in the entirely rearranged electronic cluster core. In Fig. [3] we compare 
the obtained dependence f3{uj) for the d-shell photoionization of the Na^g cluster. As in 
the case of the Na^ cluster different experimental dependencies B,C,D,E correspond to 



8 



1,5 
1,0 
0,5 
0,0 
-0,5 
-1,0 





• 


» « >• J, ▼ ▼ 

j , . 


r-form (RPAE) 

v-form (RPAE) 

■ B[2] 

• C[2] 
D[2] 

T E[2] 



1 2 3 4 5 

photon energy, eV 

FIG. 3: Angular anisotropy parameter /3 for the photoionization cross section of the Na^g cluster 
versus energy of photon; solid and doted lines are the RPAE results; dependencies B,C,D,E are 
the experimental data 1|. 

the photoionization from the sublevels of the 3rf-orbital sphtted by the crystaUine field of 
the cluster. Note that the calculated dependence of on uj possesses a characteristic 
minimum in the vicinity of the threshold, which is well seen in the all experimental curves. 
A small deviation of the parameter /J calculated using the length and the velocity form of 
the amplitudes Aii arises due to accounting for the electronic cluster core rearrangement 
in the static approximation. A small shift in uj of the calculated curves with respect to the 
experimental ones is a result of some deviation of the experimental value of Ip from the one 
calculated within the jellium model framework. 

In this paper for the first time the consistent many-body theory has been applied to the 
calculation of the angular distributions of photoelectrons emitted the process of photoion- 
ization of metal clusters. The concrete analysis was performed for sodium cluster anions 
Nttj and Naig possessing nearly spherical geometry arising due to the closure of all elec- 
tronic shells of delocalized electrons. For these clusters we have calculated the parameter of 
angular anisotropy /3 of the photoionization cross sections. These calculations revealed the 
very important role of many-electron correlations in the formation of angular distributions of 
photoelectrons and explained the behavior of the angular anisotropy parameter l3{uj) versus 
photon energy for sodium clusters recently measured in 

The next logical step of these investigations could be an extension of the formalism pre- 
sented in this Letter for clusters with partially filled shells. This treatment is to be performed 
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within the framework of the deformed jelhum model llj. The most detail description of 
the problem, of course, should take into account the detail geometrical structure of the ionic 
core. Also, for the accurate description of the near threshold behavior of the photoionization 
cross sections of negative ions one need to account for the effect of dynamical polarization 
of the electronic core, i.e. for the polarization of the electronic core by the outgoing elec- 
tron, as well as the relaxation processes of the electronic core. This can be achieved by the 
method combining the Dyson equation method and RPAE [g], [l^ . All these problems can 
be addressed for different type of metal clusters (Au, Ag, Mg, Al, K, Sr etc) and remain 
open for further investigations in this field. 

We are thankful to Dr. Andrei Korol and Dr. Andrei Ipatov for helpful discussions. 
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